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PACS 05 . 45 . -a - Nonlinear dynamics and chaos 

PACS 71.55. Jv - Disordered structures; amorphous and glassy solids 

PACS 42.65.Sf - Dynamics of nonlinear optical systems; optical instabilities, optical chaos and 

complexity, and optical spatio-temporal dynamics 
PACS 37 . 10 . Jk - Atoms in optical lattices 

Abstract - Bose-Einstein condensates loaded into kicked optical lattices can be treated as quan- 
tum kicked rotor systems. Noninteracting rotors show dynamical localization in momentum space. 
The experimentally tunable condensate interaction is included in a qualitative Gross-Pitaevskii 
type model based on two-body interactions. We observe strong and weak chaos regimes of wave 
packet spreading in momentum space. In the intermediate strong chaos regime the condensate 
energy grows as t^^'^ . In the asymptotic weak chaos case the growth crosses over into a t^^^ law. 
The results do not depend on the details of the kicking. 
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Introduction. — There is much interest in quantum 
behavior of systems that are chaotic in the classical limit. 
The hallmark model used in many such studies is the 
kicked rotor [l][2], which was observed to exhibit a para- 
metric transition from integrable dynamical behavior to 
a fully chaotic one. In physical terms, this corresponds 
to an unbounded diffusive increase of the rotor energy in 
the chaotic regime, as opposed to energies bounded by in- 
variant tori in the integrable case. Once quantized, the 
kicked rotor disrespects the classical chaotic energy in- 
crease - a highly interesting phenomenon known as dy- 
namical localization, similar to Anderson localization in 
disordered tight-binding models [sjjl]. Dynamical local- 
ization has been observed in laboratory experiments with 
atoms in microwave fields i5j|6j, and in ultracold dilute 
gases condensed within kicked optical lattices [7]-[9]. In 
particular, the latter realization with a Bose-Einstein con- 
densate, where the magnitude and the sign of the interac- 
tions between atoms can be tuned via Feshbach resonance 
by an additional magnetic field, has greatly opened op- 
portunities to study dynamical systems in the presence 
of many-body interactions. In the mean-field approxi- 
mation, these interactions can be modeled by adding a 
nonlinear cubic term in the Schrodinger equation, giving 
the well-known Gross-Pitaevskii equation. This can also 
be extended to the quantum kicked rotor model, result- 



ing in its nonlinear generalization. Considerations of non- 
linearity's effects on dynamical localization have initially 
shown that with strong enough nonlinearity, quantum lo- 
calization is destroyed and subdiffusive spreading occurs 
similar to results from nonlinear disordered me- 
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dia. The possibility to view dynamical localization as a 
counterpart to Anderson localization encourages further 
detailed numerical studies of this problem - particular mo- 
tivation is found in recent studies ,12,^13^ whereby a quan- 
tum kicked rotor multiply kicked with three incommensu- 
rate frequencies mimics the three-dimensional Anderson 
metal-insulator transition. In this paper, we shall focus 
specifically on the application of a theoretical framework 
developed for spreading of waves in nonlinear disordered 
media jT4 17 . Our main aim is to establish spreading 
laws for a nonlinear version of the quantum kicked rotor. 
In the first section, we give the model description. Next, 
theoretical predictions of spreading regimes are discussed. 
Numerical results are then presented, followed by our con- 
clusions. 

Model. — Experiments of quantum kicked rotor sys- 
tems within Bose-Einstein condensates, where many-body 
interactions play a significant role, focus theoretical atten- 
tion on dynamical localization in the presence of nonlinear 
interactions. In the mean-field approximation, dynamics 
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of the kicked rotor can be modeled by the following form 
of the Gross-Pitaevskii equation 

= - 1^ 2 + /3 I^N + fc cos(^) • V E - ^^)- 

(1) 

Here {3 is the nonlinear strength, which is proportional 
to the tunable two-body scattering length of atoms in a 
BEG. M is the mass of the atoms, k is the perturbative 
kick strength, and T is the period of applied kicks. Note 
that the analogy between an abstract rotor and the atomic 
wavefunctions is obtained when the atoms are loaded 
into a momentum eigenstate of the lattice with a Bloch 
wavenumber of zero - spatially homogeneous kicks will 
keep the Bloch wavenumber invariant, yet allow changes 
in the momentum (see below). 

When the nonlinear strength is zero (/? = 0), the prob- 
lem is reduced to a linear map. Because of the instant 
action of the perturbation, the evolution operator U of 
the linear system can be written as product of two non- 
commuting unitary operators 

^{6, t + l)^U ■ tl){e, t) = i3(fc, 6) ■ G(r/2, 9) ■ ij{0, t) 

The first corresponds to free rotation between two succes- 
sive kicks, 

G(r/2,.)=exp . . -, 

and the second describes the evolution over the kick, 

B{k,d) = cxp{—ik cos 9) , k = k/h. 

The solution ip{9, t) can be expanded in an angular mo- 
mentum basis 

_j oo 

^(^'*) = ^ E (2) 

^ n— — oo 

where the coefficients An{t) are Fourier coefficients of the 
time-dependent wave function ip{9,t). As a result of the 
action of the evolution operator U on 'ip(9,t) over one pe- 
riod T, the following map for the Fourier coefficients is 
obtained 

A„{t + 1) = J2i-'T'"''Jn-ra{k)Ara{t)e~'^"'' (3) 
m 

where J„_m(fc) is a Bessel function of the first kind. Un- 
like the classical rotor, the behavior of the quantum one 
essentially depends on two parameters. The perturbation 
strength k sets the effective number of unperturbed states 
covered by one kick, while r relates the period of applied 
kicks T to the natural frequency of the rotor, which is de- 
fined a,s u! — h/2M. By w = 2tt/T, the parameter t can 
be written as t = AnT/f. 

Since the perturbation is time-periodic, Floquet theory 
can be applied flSl 

i,{9,t) = e-'^'^^i9,t), M(^,t + 1) = 0^(0, t) (4) 



where 4>x{(^i ^'"^ t-periodic functions with period T (Flo- 
quet states), and x are the quasienergies. Then Eq.([3]) can 
be treated as an eigenvalue problem 

X.A^^ = J2i-^r-"'Jn-m{k)e-^^"''A^^ (5) 

m 

with eigenvalues located on the unit circle in the com- 
plex plane: X^, = e'^". For a rational ratio of T/T, ex- 
tended eigenvectors are obtained. This case corresponds 
to the regime of quantum resonances, where the rotor en- 
ergy E{t) = J2n\^n\'^ /"^ grows as t'^ (balfistically) in 
time - the cause being quantum interference effects 
For an irrational ratio of T/T (the case treated here), the 
sequence {e~*5'" \ m G Z} is quasiperiodic and the cor- 
responding eigenvectors are exponentially localized - this 
effect is dynamical localization. If the argument Tm^/2 
in the exponential is replaced by a truly uniform random 
sequence in the interval [0, 27r], calculations have shown 
the same effect of exponential localization in correspond- 
ing eigenvectors, supporting the idea that quasiperiodic 
sequences mimic disorder [19] . 

For a given set of parameters r and k, the localization 
extent of an eigenvector is estimated by its participation 
number: P = l/X^rJ^nl^- define the localization 

volume as the average of this over the mode index, V — 
{P)y- This localization volume depends on the strength 
of kick k, and is increasing with k. It measures the typical 
size of a localized eigenvector. The parameter r has no 
significant influence on the localization volume, provided 
that the ratio T/47r is irrational. The localization volume 
also deflnes the number of other modes a given eigenvector 
can interact with, once the nonlinearity is added. Since 
all these V modes have in general different quasienergies 
Xjy, we can introduce a characteristic difference scale - the 
average quasienergy spacing d. Since all quasienergies are 
modulo 27r, it is simply given by 

d = 2t:/V . (6) 

In the nonlinear case (/? ^ 0), expanding ip{9,t) in the 
angular momentum basis, Eq.([2]), the dynamics between 
two successive kicks is described by 

= ^2^""^^" ^ ^E E ^rii^n2^T!-(n2-ni)' C'') 
rii n2 

where /3 — f3T/2TTh. Since we are ultimately interested 
in the limit of small amplitudes, we treat the nonlinear 
term perturbatively and keep only the diagonal terms in 
Eq.([7]). Neglecting all off-diagonal terms in Eq.([7| and 
integrating over the free motion between two delta kicks, 
An{t) evolves according to 

A„(t+l) = A„(t)e-'i"'+^'5|^"l', (8) 

After additional integration over the infinitesimal interval 
of one kick, the map - which now describes the evolution 
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over one whole period - becomes 



with the overlap integral 



Anit 



'l/,/J-l,^2,M3 



(9) 

This map was first introduced by Shepelyansky in [To] . 
Comparison of the results of this map with direct numer- 
ical simulation of the corresponding model, Eq. ([!]), has 
shown differences on a short time scale, but the same 
asymptotic behavior in the rotor energy 20 . At the same 
time, this model allows for more efficient and faster nu- 
merical computation. 

Numerical simulations of the map, Eq.([9|, are per- 
formed for different parameter sets {fc, r, /?}, starting with 
a single-site excitation at site no: = S^ na- The 

lattice sizes were between N = 2^° for smaller k values 
and N = 2^^ for larger values. To characterize the time 
evolution properties, we compute the on-site probability 



(12) 



Since all eigenvectors are exponentially localized, each is 
effectively coupled with a finite number of neighbors - the 
interaction has finite range. However, according to norm 
conservation, as the wavepacket spreads, individual on-site 
probabilities will start to reduce - so will the average norm 
density, n. Since the coupling strength between eigen- 
modes is proportional to n, it also decreases with packet 
spreading, while the number of excited eigenmodes grows. 
In disordered nonlinear wave equations, spreading 



regimes were observed 16 17 to depend on relations be- 



the time-dependent participation numbers 



w„, and the average energy. Since localization ^he wavepacket is self-trapped (2T1I22 



Wn = \A„ 

is in momentum space, the average energy is defined by 
the second moment of the on-site probability, multiplied 
with one half 



tween a nonlinear- induced frequency shift, the linear spec- 
tral width, and the average frequency spacing between 
interacting eigenmodes. In the case of frequency shift 
larger than the linear spectral width, at least a part of 

In weaker cases. 



(10) 



where n — nwn is the first moment. To quantify the 
sparsity of a packet we calculate the compactness index as 
a ratio of the square of the participation number and the 



second moment, C 
ics, we ensure the total norm 
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absolute relative errors of < 10" 



. To validate our numer- 
= Wn, is conserved (at 
^) during all simulations. 



Regimes of Spreading. — In Eq.Q, the nonlinear- 
ity induces a quasienergy shift of A0,„ = It 
also induces an interaction between localized eigenvec- 
tors. Namely, the nonlinear term of the map in Eq.([9| can 
be rewritten as exp(i/3 |^,„|^) = 1-1- /(«/? |A„i(t)|^). Ex- 
panding An in the linear eigenbasis, An{t) — 4',^{t)A^, 
Eq.Q is transformed to 

Mt + 1) = A,0,(t) + ^ Mt) ■ {-lY-^Jn-m{k) ■ 



This is the corresponding map in the basis of localized 
complex eigenmodes, (j)y{t). Under small quasienergy 
shifts, the nonlinear term is Taylor expanded (/ ~ 
il3\Am{t)\^) to yield 



<^,(i + l) = A,0,(t)+/3 ^' 



'P^n(f>*a2'^P■3 (11) 



,AI2 ,M3 V^A'l t'Ms 



when self-trapping is avoided, two outcomes are possible. 
If the frequency shift is larger than the average spacing, 
almost all eigenmodes resonantly interact with each other 
- the strong chaos regime. Otherwise if the frequency shift 
is less than the average spacing, there is very minimal res- 
onant interaction amongst eigenmodes - the weak chaos 
regime. 

The similarity of the models hints that the spreading 
conjecture of 14 may also be valid for the nonlinear quan- 
tum kicked rotor. In this conjecture, an exterior eigen- 
mode within a localization volume of the wavepacket 's 
boundary is incoherently heated by the packet. The 
probability of a resonance inside the packet depends on 
the ratio between the nonlinear-induced quasienergy shift, 
S « /3n, and the average quasienergy spacing d: V{Pn) « 

~l3n/d 



1 _ g-P"/'' [14)[15]. From Eq.([llj), it follows that the ex- 
terior eigenmode is excited according to 

Mt + 1) « K(^At) + p7v"^v{MIii), (13) 

where f{t) is a stochastic force. The time-dependent ex- 
terior eigenmode norm is \(t>{t)\ ^ l3'^-n?{V{l3n)yt. It ex- 
cites to the level of the packet norm density n in a time 
T ^ l3-^n-^{V{Pn)y^. At this point, the once exte- 
rior eigenmode can be considered to be engulfed into the 
packet interior, so the rate of norm diffusion can be de- 
fined as D = T^^ ^ l3'^'n?{V{l3n)Y . The second moment 
is inversely proportional to the square norm, m2 ~ 1/rt^ 
and the diffusion equation is m2 Dt. From these two re- 



lations, one obtains an equation \/n 



(1- 



Given the limits of the resonance probability, this yields 
two regimes of subdiffusive spreading 

m2 ^ Pt^^^, (in/d > 1, strong chaos 

m2 - ^"2/3/3^/3*1/3, I3n/d<l, weak chaos 

The average energy of the quantum kicked rotor E is 
defined by the second moment of the on-site probabilities. 
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Eq.(lO). Hence, from the two above regimes, theoretical 
expectations for the subdiffusive growth of the average en- 
ergy, E ^ t", are a « 1/2 in the strong chaos regime, and 
a ~ 1/3 in the weak chaos regime. Spreading increases 
the wavepacket width, but under a norm conservation, 
individual on-site probabilities are reduced and the non- 
linear quasienergy shift likewise decreases. A wavepacket 
launched in the strong chaos regime consequently spreads 
to a point whereby quasienergy shifts are reduced to only 
minimal resonant interactions. The strong chaos regime 
is therefore a transient one; eventually the dynamics will 
cross over to approach the weak chaos regime in an asymp- 
totic limit. Nonetheless, the previous studies in disordered 
nonlinear Anderson or Klein- Gordon chain have shown 
that this transient regime can exist for many orders of 



magnitude before the crossover is observed 16 17 



Numerical Results. — Since all eigenvalues are lo- 
cated on the unit circle for the quantum kicked rotor, there 
is essentially no limit to the linear spectral width. There- 
fore, unlike in disordered chains, the self-trapping regime 
is not expected to be observed, regardless of the nonlinear- 
ity strength. For irrational r, all quasienergies are basi- 
cally distributed homogeneously over the interval (0, 27r). 
The localization volume, V, depends on the kick strength 
k. The average spacing between quasienergies can then be 
estimated using Eq.(|6]). In Figjl] the dependence of the 
average spacing d on the kick strength k is shown, for a 
quasiperiodic sequence (t = 1). The dependence d{k) for 
a random sequence is quite close (well inside the shown 
standard deviation) to that obtained for the quasiperiodic 
case, therefore it is not shown. 



10 



10 



to 

o 
a 

3 

cr 



10' 



o 
Z 



10' 



: 1 


1 

+ 


II: 

Strong Chaos : 


Weak Chaos 

1 


+ 

1 


+ 

1 , 1 



3 4 
Kick Strength, k 



Fig. 1: The parameter space as given by the nonlinear 
quasienergy shift S versus the kick strength k. The sepera- 
tion between the two expected regimes is estimated by d, the 
average spacing between quasienergies. The full line is d for 
the quasiperiodic sequencing (r = 1), while dashed lines cor- 
respond to its error within one standard deviation. Crosshairs 
correspond to parameters used for numerical simulations. 



Starting from a single-site excitation, before quantum 
supression effects become significant, the packet will ex- 
tend over the localization volume. Consequently, its norm 
density will drop from the initial value Uin = 1 to n w 
1/V. Choosing values of the nonlinear coefficient /? gives 
quasienergy shifts of (5 ~ either above or below the 
line d{k). Considering that the shift 6 and the average 
quasienergy spacing d are both inversely proportional to 
the localization volume, the condition 5 > d [k) can be 
reduced to /3 > 27r. This effectively gives us a param- 
eter space in which to choose the nonlinearity parame- 
ter p. The weak chaos regime is expected for the values 
/3 < 27r, while the strong chaos regime is expected for 
values /3 > 2tt. Simulations were performed for the val- 
ues fc = 3, 5 for both a quasiperiodic sequence (r = 1) 
and for a random sequence. For each set of parameters, 
the simulations were over more than 500 different realiza- 
tions. For random sequences, these realizations were sim- 
ply unique uncorrelated sequences, while for the quasiperi- 
odic case the unique realizations were different initially 
excited states of the quantum kicked rotor. To suppress 
fluctuations, the logarithms (base 10) of time dependent 
measures are averaged over these realizations - we do this 
for the average energy in Eq.(lO) to give (logj^pi?). In 
order to quantify our findings, we utilize locally weighted 
regression smoothing ^3j, and then calculate the corre- 
sponding local derivative 



a 



d{log,oE) 



(14) 



d (logio t) 

via a central finite difference [24]. 

For fc = 3, the results are seen in Fig[2j In the weak 
chaos regime, /3 = 0.4 (blue curves), the exponent a in- 
creases slightly above the expected 1/3, but reaches it 
at later times. In the strong chaos regime, (3 = 12 (red 
curves), the exponent a decreases from the predicted value 
1/2, without spending much time in the strong chaos 
regime. These findings are then compared with the re- 
sults for random sequences, shown as dashed lines. As 
can be seen, spreading behaviors are quite similar to 
the quasiperiodic sequencing. Additionally in the inset 
of Fig[2) the average compactness index (C) asymptoti- 
cally attains values in the range ^ 3.4 — 3.6 for both the 
quasiperiodic and random cases. Note that the compact- 
ness index does not significantly decay, meaning that sub- 
diffusive growth of the second moment is followed by cor- 
responding subdiffusive growth of the participation num- 
ber. This indicates almost all states inside the packet are 
homogeneously populated during the spreading 11 . 

Lower values of the kick strength k correspond to 
smaller localization volumes. These lower values may be 
associated with no clear identification of the strong chaos 
regime in the kicked rotor; since small localization volumes 
are related to a small interaction range between eigen- 
modes, it ought be that the strong chaos regime is quite 
short in time for lower values of the kick strength k. To 
test this argument, we performed simulations for k = 5, 
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Fig. 2: Under a kick strength of = 3, measures for /3 = 0.4 
(blue) and 13 — 12 (red), for both quasiperiodic sequences set 
by r = 1 (solid line), and for random sequences (dashed line). 
Upper row: Mean logarithms for energy (log-^gE). The clouds 
around the quasiperiodic sequences correspond to one standard 
deviation error. Lower row: finite-difference derivative of the 
above. Grey horizontal lines correspond to exponents for weak 
and strong chaos regimes. Inset: average compactness index 
{Q as a function of time. 



a stays for a significantly long time near its theoretically 
predicted value of 1/2. The compactness index shows the 
same behavior as in A; = 3, indicating homogeneity of the 
populated states within spreading packets. The results for 
the random sequences (dashed lines) also confirm similar 
behaviors for k — 5. In particular, at the end of the ob- 
served time, it can be seen that the exponent a decreases 
below 1/2, signaling an entry of the system into a regime 
crossover from strong to weak chaos. 




log t 
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but with the total number of kicks one order smaller than 
for fc = 3. The reason for this reduction is simply that 
bigger k requires larger system sizes, increasing computa- 
tional time. 

For fc = 5, the results are shown in Figj3j The deriva- 
tive a, for /3 = 0.3 (blue curves), shows similar behavior 
to the case for fc = 3 in the weak chaos regime - namely, 
it grows until it slightly exceeds and then drops back to 
1/3. Even though for fc = 5 we are unable to look at 
extremely long times, we can still conclude that observed 
behavior is not much qualitatively different in the weak 
chaos regime than for the case fc = 3. On the other hand, 
for /3 = 10 (red curves), which should correspond to the 
transient strong chaos regime, it is clear that the exponent 



Fig. 3: Under a kick strength of fc = 5, measures for /3 = 0.3 
(blue) and P — 10 (red), for both quasiperiodic sequences set 
by r = 1 (solid line), and for random sequences (dashed line). 
Upper row: Mean logarithms for energy (logj^Q E) . The clouds 
around the quasiperiodic sequences correspond to one standard 
deviation error. Lower row: finite-difference derivative of the 
above. Grey horizontal lines correspond to exponents for weak 
and strong chaos regimes. Inset: average compactness index 
{() as a function of time. 



Conclusion. — The results of extensive numerical 
simulations show that in the presence of nonlinearity, dy- 
namical localization is destroyed. Techniques of averag- 
ing and local derivatives on logarithmic scales indicates 
existence of the theoretically predicted spreading regimes 



p-5 



G. Gligoricl,2 J. D. Bodyfeltl S. Flachl 



of weak and strong chaos, as well as the crossover be- 
tween them. This supports the idea that delocalization 
can be understood as a result of mode interactions inside 
the packet with those outside, as induced by the nonlin- 
earity shift. In addition, the existence of a crossover might 
be associated with findings that the exponent of subdif- 
fusive growth of the second moment can be temporarily 
close to 0.4, which was obtained in earlier papers on the 
nonlinear quantum kicked rotor 



10 ,11. Our simulations 



also confirm a novel lack of a self-trapping regime, which is 
not expected in this model due to an unbounded - yet cir- 
cular - linear spectrum. We anticipate higher-dimensional 
generalized theoretical frameworks developed for nonlin- 



ear disordered chains 14 to also be valid in the "higher- 
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